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Preceding  page  blank 


ON  SOME  DISTRIBUTION  PROBLEMS  CONCERNING  THE 


CHARACTERISTIC  ROOTS  OF  S^"1  UNDER  VIOLATIONS 


by 

K.C.S.  Pillai  and  Sudjana 


,1.  Introduction.  The  distribution  of  the  characteristic  roots  of 

under  certain  violations  was  obtained  by  Pillai  [12].  In  the  present  paper 

Pillai' s  distribution  has  been  used  to  derive  the  following:  1)  the  density 

2 

function  of  T,  a  constant  times  Hotelling's  TQ,  2)  the  moments  of  T, 

3)  the  moment  generating  function  of  T,  4)  the  m.g.f.  of  Pillai's  trace, 

5)  the  distribution  of  Wilks'  criterion  and  6)  the  density  function  of  the 
largest  root  in  two  forms.  These  results  are  useful  in  studying  the  exact 
robustness  of  at  least  two  multivariate  hypotheses,  namely,  a)  MANOVA  under 
violation  of  the  assumption  of  a  common  covariance  matrix  and  b)  equality  of 
covariance  matrices  of  two  p-variate  normal  populations  when  the  normality 
assumption  is  disturbed.  The  exact  robustness  comparisons  will  be  reported 
later. 

2.  The  density  function  of  T.  In  this  section  we  will  derive  the  density 
function  of  T  =  X  tr  where  (p  x  p)  is  distributed  WCpjn^Lj.fl) 

i.e.  non-central  Wishart  distribution  on  n^  d.f.  with  non-centrality  Q  and 
covariance  matrix  and  S2  (p  x  p)  independently  distributed  central 

Wishart  W(p,n.,Z„,0)  where  n.  ,n-  >_  p  and  X  is  a  positive  real  number. 


In  the  case  of  n,  <  p  w.  have  only  non-zero  roots  of  S^1  and  the 

density  function  of  T  can  be  obtained  from  that  for  n,  >  p  if  ™  the 
latter  case  the  substitutions  (nj.nj.p)  +  CP -n i*n2"p ,n2^  are  “de' 

Let  us  write  A,  •  |  I,  2  ^  I,  2  and  *2  *  T  !l  *  -2  ’  the" 

applying  (29)  of  James  [6]  we  get  the  joint  density  of  A,  and  A,  as: 


Cl(P.n2,  e  -l^l2  2  •  - 


tr  Z  -  -5^1 
e  -  |Z|  2  1 


Rc(Z)-X0 


exp(-tr(I-W)A^)  dZ 


ip(p-i)  ^>(P*D  i 

where  CjCp.nj)  -  2  /[(2»i)  Pp  2 


I'p(f’2)1  ’ 


A  -  t\  S'1  t\  .  W  •  02  2'1  »2  »«•  l  *  *0  *  1  l  ”lth  -0  P’d>  Sym”etriC 
Mtril  and  Y  a  non-singular  real  symmetric  matrix  such  that  (X-»)  non‘ 
singular.  Not.  that  the  roots  are  invariant  under  the  above  transformations 

7  2 

and  also  under  the  following  transformation  Bj  -  (I-W)  ^(I-W) 

1_  1 

B  -  (I-W)2  A2(I-W)2.  Using  these,  we  obtain  the  joint  density  of  ^  and 
*2  “*  **  ** 

Bj  as  follows: 

■tr  A  in2  -tr?l|„  ,7("rP'I)|B  i(Vp'1) 

(2.1)  C1(p,n2)  e  |a|  e  '-2* 


tr  Z  -  ’  2^ni+n2^ 

.  /  e  -|Z|  2  l|l-?l 


_  1  .  1 
expC-tr(I-W)  2A(I-«) 


Re(Z)-X0 


Now  the  Laplace  transform  of  T  ■  X  tr  =  X  tr  is  given  by 

E(exp(-t  X  tr  B^-1)).  So,  to  find  this  transform  multiply  (2.1)  by 


2 


,-l, 


expC-t  X  tr  )  and  integrate  Bj  out,  we  get  the  Laplace  transform  of 
T  in  the  form: 


r,  ..  -tr?,.|J"2/  2»*1 

C2(p,n1,n2,X)e  |A|  t 


/  / 


tr  Z  -bn  -  i(n.+n  )  i(n  ♦n.-p-l) 
e  ~  |  Z  |  2  Ml-Wl  2  1  2  |B  |2  1  2 


B2>0  Re(Z)-X0 

1^  .1  .1 

•  1 1  ♦  (tX)_1B2|  2  1  exp(-tr(I-W)  2A(I-W)  2B2)  dB2  dZ  , 


^P(P-I)  i  4p(P+l)  , 

where  C2(p,nlfn2lX)  -  2Z  r  C^)/ [(2iri)z  X*  rp(±n2)]  . 

1  1  II 

T  2  "22 

Now  letting  D2  ■  A  (I-W)  B2(I-W)  A  ,  the  above  expression  becomes: 


-tr  0  -  ■=«. 

(2.2)  C  (p.n  n  x)  e  '|a|  /  e 

d2>o 


_tr  °2  ,2(ni+n2"P-1) 

l?2l 


Re(Z)-X0 


_  1  1 
tr  Z  -  i  i 

e  "|Z|  '*  1 1 1  +  (tXj'^I-WJA-1  D2|  1  dD2  dZ 


Further,  let  us  write 


-  b\ 

1 1  ♦  (tX)-1  (I-W)A  '1D2|  2  1 


;  ,  4  -  .  .,-k  v1'1?-' !)*■*  y 
J0  l  ’  in 


in  (2.2)  and  integrate  term  by  term  with  respect  to  t  for  sufficiently 
large  values  of  Re(t)  and  then  integrate  out  D2  making  use  of  (12)  of 
Constantine  [1)  to  obtain  the  density  of  T  in  the  form: 


3 


-tr  ft 


■hi 

^(p.n^nj.XJT  (i(n1+n2))e  ~ | A |  2  1  T 


-1 


/ 


trZ,  r  r  (Ini^4(ni+n2»lc(-T)k  ] 

e  •'|Zj  2  1  I  I  --1"  y-«  2  - -  C  [(XA)_1(I-W)]dZ 

Re  Zs,x^  -  ’ -  "  * 


k»0  k  kl  r  (^piij+k) 


which  is  convergent  for  |  T/XX ^  J  <  1,  where  X^  is  the  minimum  root  of  A 
(see  (8]).  Now  (XA)"1  being  symmetric,  can  be  diagonalized  by  an  orthogonal 
transformation  H  e  0(p)  and  perform  the  integration  over  0(p),  then 
finally  apply  (17)  of  Constantine  [2]  to  get  the  density  function  of  T  as: 


i  i  -  ^n,  -tr  ft  -rpn.-l 

(2.3)  f(T)  *  [rp(|-(n1+n2))/rp(in2)]|XA|  2  1  e  ~  T2  1 

1  k  -1  -1  T(ni  ”P~1) 

.  i  i  eg) 

k«0  K  kl  Tptjfryk)  C^(I) 

Formula  (2.3)  will  give  special  cases  as  follows: 

a) .  For  ft  -  0  which  implies  (0)  =  ( Jn  )  C  (I)  where  y  =  I(ni-p-l), 
we  have  the  result  of  Khatri  [8]  formula  (9) . 

b)  .  For  A  *  I  and  X  =  1  we  have  Theorem  4  of  Constantine  [2], 


3.  The  moments  of  T.  In  order  to  obtain  the  moments  of  T,  we  note  that 

T  3  (X  tr  I  Sc(BlB2^*  w*iere  Bi  and  B2  are  as  Section  2 

and  we  shall  use  the  joint  density  of  B:  and  B2  in  (2.1).  The  kth  moment 

of  T  is  E(T  ).  Multiplying  (2.1)  by  X^  £  and  integrating  out 

< 

first  B2  using  (8)  of  Khatri  [7]  and  then  Bj  using  (12)  of  Constantine  [1] 
we  have: 


i 
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E(Tk)  =  C1(p,n2)  Xk  e 

.  /  '  |  Z  | "  ^U-Wf  l  r  (ln2r<)r '  (4b  ..k)C  [Ad-IO'^dZ 

Re  Z=XQ  '  ~  ~  k  F  F 

Now  transform  Jl  +  HAH'  where  H  e  O(p)  and  integrate  over  0(p),  then 
the  right  hand  side  in  the  above  expression  becomes 

-tr  fl  j  j  CJXA) 

I  rpf2ni»K^rp^2n2,"'C^  C  (1^ 

K  K 

tr  Z  —  -  ^ i 

.  /  e  **  |  Z  |  2  1 1 1  —  W |  2  1  C  f(I-W)  X)  dZ 

Re  Z*XQ 

I  i 

Replacing  W  by  ft2  Z_1  fl2  and  noting  that  Ck((I-W)  *)  =  (I  +  n(Z-£l)  *) 
and  rearranging  the  necessary  variables, we  obtain 

v  l  1  CJXA) 

E(Tk)  =  c1Cp,n2)  l  r  -c-rry 

K  “  r  K\ 

1 

tr  U  -  *n.  , 

•  /  e  ~ |U|  Z  1  C  (I+nU  A)dU  , 

Re  U>0  ~  K  ~ 

where  U  =  Z  -  n  Now  using  (17)  of  Constantine  [2]  we  have  the  kth  moment 
of  T 

-(n  -  -1) 

(3.1)  E(Tk)  =  [rpC^2))_1  I  rp(|n2,-K)CK(XA)L2  ">  P  (a)/cK(i)  . 

Finally  using  the  fact  that  r  (t,-<)  =  [(-ljr^(t)]/(-t*  j(ptl))K  ,  we  get 
kt^1  moment  of  T  as 
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(3.2) 


» 


nk,  ,  ^k  r  CK™  L< 


yCnj-p-l) 


(-ft) 


E(TK)  =  (-1)K  I  4 

<  4(P+l-n2))K  Ck(I) 

which  exists  only  for  n 2  >  2k+p-l. 

Formula  (3.2)  will  give  special  cases  for  special  values  of  ft  and  A. 
If  we  let  fl  ■  0,  (3.2)  gives 


E(Tk)  =  (-l)k  l  (^n1)K  Cic(AA)/(|(p+l-n2))ic  , 


which  is  the  result  of  Khatri  [8] ;  except  that  his  formula  contains  an  error 
in  the  denominator.  His  denominator  is  (y(p-n2-l))j  and  the  right  one  is 
(■j(p-n2+l))j.  Substitution  of  A  *  I,  A  *  1  in  (3.2)  will  give  the  result  of 
Constantine  [2]  formula  (38) . 

4.  Moment  generating  function  of  T.  Pillai  [12]  has  obtained  the  joint 
density  function  of  the  roots  r. , . . . ,  r  of  S.S”1  which  has  the  form: 

i  p  >1.2 

-tr  fl  -  y(n.-p-l) 

(4.1)  C(p,n  ,n  )e  ~ | A |  *L\r\Z  n(r.-r.) 

A  -  -  i>j  3 


c  (ft) 


I  I— 

k=0  Kfinp^^pk!  n=0  v;s 


;  r  t-1’" 

i  1  - - C”0)  n! -  » 


where  g  are  constants  (Constantine  [2],  Pillai  and  Sugiyama  [14]), 

<,  v 

R  »  diag  (Tj,...,  Tp)  with  0  <  Tj  <  ...  <  r  <  »  ,  +  ...  +  6p  *  k  +  n, 

6  «»  (6^, . . .  ,<5p)  and 

1  2 

(4.2)  C(p,nrn2)  =  ^(n^))/ [  rp  (^nA )  Tp  (|n2)  Tp  (^p)  ]  . 
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To  obtain  the  m.g.f.  of  T  ■  Atr  SjS^1  we  shall  take  the  expected  value  of 
exp  (t  A  tr  R)  with  respect  to  (4.1),  since  under  the  transformations  employed 
by  Pillai  [12]  to  get  (4.1),  tr  SjS’1  ■  tr  R.  Now  transform  back  ft  +  H  R  H' 
where  H  e  O(p)  and  the  latter  matrix  R  is  symmetric,  then  perform 
the  necessary  integrations  (i.e.  use  (44)  of  Constantine  [1],  and  integrate 
out  R  >  0) ,  we  get : 


-tr  ft  -  4n. 

E(exp(t  tr  R))  -  C^p, n^n^e  ~|a| 


”  I  i  (-i}  avvjd 

k-0  *c  k!  Ck(I)  n-0  v,6  nI 


*  / 

R>0 


tr(tAR) 

e 


|R|  1  c6(A  *r) 


dR 


where  Cjfp.nj.n^  »  rp Cy(n1*n2) ) / [ Tp (^sij )  Fp (^ri2)  j 

Let  us  now  apply  (12)  of  Constantine  [1].  IJpon  simplification  we  finally 
obtain  the  m.g.f.  of  T  as: 

_  1 

1  1  "tr  n  ” 

(4.3)  E(exp(t  A  tr  R))  *  [rp(J(n1+n2))/rp(ifi2)]  e  |-t  A  A| 


l  l 


ckW 


>n,l 


k-0  k  (in.)  k!  C  (I) 
2  I  K  K  _ 


1  <1 


-(tA) 


■V1) 


n»0 


Again,  as  before  we  can  put  special  values  of  ft  and  A  to  get  special 
cases.  Making  a  substitution  ft  *  0  in  (4.3)  and  noting  that  g^  v  =  1, 

6  =  v  we  have: 


- 

(4.4)  E(exp(t  Atr  R) )  -  Crp (^■fn1*n2) ) /rp(^n2) ]  | -t  A  A|  2  1 

2F0(I<V  ^ni+n2);  (tx)_1 
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while  letting  A  -  I  and  X  ■  1  and  v  ■  0  which  implies  <  ■  d  we 

obtain: 


1 


1, 


(4.5)  E(exp(t  tr  R))  -  [r  (~(n^n2))/r  (^J ]  (-t) 


*  -tr  G 

L  M  /•  *  \  ‘  A  _ 


•iF0(j(ni*n2);  -t'1  n) 


or  in  an  alternative  form: 


‘  I*"), 7(VV 


(4.6)  E(exp(t  tr  R))  »  [T  (•|-(n1*n2))/r  (-|«2)  ]  (-t)  z  J 1 14-t" An 


5.  Moment  generating  function  of  V^.  To  obtain  the  m.g.f.  of  Pillai's 
criterion  which  is  defined  by  *  tr[(XR) (I* XR)"^]  where 

R  ■  disg  (Tj,...,  Tp),  Tj,...,  Tp  being  the  roots  of  SjS"1,  we  start  from 
the  following  joint  density  of  roots  of  SjS"1  (Pillai  [12]): 

-tr  fi  -  ii.  4(n. -p-1)  -  -kn.+n  ) 

C1(p,n1,n2)e  -|a|  ^Rl2  |I*AR|  n  (r.-r.) 

’  '  '  3 

tr  z  -  1|  1  1 

•  /  e  ~  |  Z|  7l  F  (i(n  ♦n  );  1-x'1  (I-W)TA-1  (I-W)7.XR(I+XR)_1)dZ 

Re  Z>0  *  i  v  £  i  £  . 


where 


CjCp.nj.n2) 


121,,.  1  , 

^9  ^P(P-l)  ,  2?(P+i)  ,  , 

"  2  r  4oij*n2))/[(2wi) 


Let  L  *  AR(I+AR)  then  the  above  density  gives  the  joint  density  of 

roots  of  L  as: 

1  P 

'tr  T^i-P'1)  y(n.,»p-l) 

Cj(p,nj,n2)e  ~  j  XA  |  A|l|2  jl-L|2  ‘  n  (£..-£.) 

"  ~  '  i>j  3 


tr  z  -  In  I  I 

/  e  ~ | Z|  2  1  F  (i(n.*n-);  I-A_1(I-W)2  A_1(I-W)2,  L)  dZ 

Re  Z>0  '  1  v  t  1  £ 


-tr  fl 
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where  0  <  l  <£„<...<  I  <1. 

1  2  p 

The  expected  value  of  exp(t.  V^)  with  respect  to  this  density  is 


-tr  n  - 

C1(p,n1,n2)  e  ~|XA| 


I  t  tr  L  I(n  -p-1)  j(n2-p-l) 

•  /  e  |lt  |i-l|  n  (l  -1  ) 

i  >n  -  -  ~  1  J 


tr  Z 

/  e  ~|Z| 

Re  Z>0 


1  ll 

2  1  lF0(7(Vn2);  I”X-1CI-W)2A"1(I-W)2,L)dL  dZ 


Now  use  (31)  of  James  [6]  and  rearrange  the  order  of  the  integration  to  get 


(d)  IPCP-U  2-P(P+1)  .  /  k 

E(exp(t  VlpJ))  -  \1T  / (2iri)2  jjXAj  2 


-  -tr  ft 

e 


-tr  S  i(n  ♦n_-p-l)  tr  Z  - 

•  I  e  -|S|7  1  2  /  e  -|Z|  2  1 

S>0  '  Re  Z>0 

I  '  I  "  I 

•  /'  exp[tr  Ctl+Sd-X-1  CI-W)2A-1(I-W)2})L] 

L>0  . 

j(n  -p-1)  j(n  -p-1) 

•  | L|  |l-Lr  dL  dZ  dS 


After  integrating  out  L  using  (47)  of  James  [6]  we  get: 


fnl  'tr  n 

E(exp(t  VIPJ))  -  C2(p,ni,n2)e  ~  |XA|  2  1  / 


tr  Z  -  •=« 

o  -Izl  2  1 


Re  Z>0 

-tr  S  j(n.+n_-p-l)  .  I  I 

7  e  *  |S|  1F1(^n1;y(n1+n2);tI*S(I-X_1(I-W)ZA'1(I-W)2))dS  dZ 


where 


For  n  ■  0,  (5.2)  givej  the  result  of  Khatri  [8]  with  a  correction  of  his 

expression  given  by  Pillai  [11].  Another  special  case  which  can  be  derived 
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from  (5.2)  is  formula  (3.5)  of  Pillai  [11]  if  in  (5.2)  we  let  A  ■  I,  X  ■  1 
and  we  use  (Constantine  [2]) 

i(n.-p-l)  .  n 

(5.3)  L‘  (0)  ■  (in  )  C  (I)  [  I 

'  s»0  o 

and  finally  we  let  6  ■  v  *  a  . 


Vi  c,w 

<KW!> 


6.  The  density  function  of  W^.  In  terms  of  the  characteristic  roots 
r^(i  ■  1,2,...,  p)  of  SjS”1,  the  Wilks'  criterion  is  defined  by 

W(P}  ■  1 1  ♦  XR|-1.  However,  in  the  derivation  of  the  density  of  W^  we 
consider  W^  ■  |l  -  L j ,  i.e.  we  let  L  ■  XR(I+XR)  1  in  the  formula  (3.7) 
Theorem  2  of  Pillai  [12].  The  joint  density  of  characteristic  roots 
lj,  *2» • • • ,  ap  of  L»  then,  is 

-tr  n  -  -L  y(n.-p-l)  y(n  -p-1) 

(6.1)  C(p,ni,n2)  e  ”  |  XA  |  l\l \*  |l-L| 


n  C^-ij)  l  l  (-4^) 

i>j  1  ■’  k-0  k 


ni+n2,  ck(!:) 


l  l  -^f 


v-x  ?  >i; 


d.0  5  (dn,),  C,(I)  Cj(I) 


where  CCp.nj^)  is  as  in  (4.2). 

To  obtain  the  density  function  of  W^,  first  we  find  the  kth  moment 
of  it  and  then  use  the  results  on  inverse  Mellin  transform  [3,4,5],  The  k**1 
moment  of  W^  with  respect  to  (6.1)  is 


11 


E(W^)h 


-tr  fi  -  Ifi, 
C(p,n1,n2)  e  ~ | XA | 

L>0 


r  r  VV  C<(L) 

I  I  (-4A  TT  CdL) 


k*-0  r. 
k 

l  l 


1 


-l.-K  ,2(nrp_1) 


a<^5  V’*  i  )  L6 


,1 


d-0  6  (|n1)<5C<s(I)  C6(I) 


n  (t.-i.) 
i>j  J 


(«) 


Now  transform  L  +  H  L  H',  where  H  e  0 (p)  and  L  is  a  symmetric  matrix. 
Integrating  out  H  using  (44)  of  Constantine  [1]  and  then  L  using  Theorem  3 
of  Constantine  [1],  we  obtain  the  h^  moment  of  in  the  form: 


-tr  n  -  ■=«, 


(6.2)  E(W(p))h  -  [rp(jOyn2))/r  (|n2)]  •  'M  2  * 


l  l 

k-0  tc 


k 

l  l 


(i'VVViVk  C^I)  i"1  r(r*bi) 


k! 


n  r(r+a.) 
i-1 


1 


a^6  C6(-X  *A  *)Lfi 


(fl) 


a 


d-Q  6  (jnl)a  Ca(I)  Ca(I) 


where  r  =  ^n2  ♦  h  -  j(p-l),  bi  c  j(i-l)  and  a.  -  j +kp _ ± ^  1  b.. 

Finally,  on  (6.2)  we  apply  results  on  inverse  Mellin  transform  [3,4,5]  to 
obtain  the  density  of  of  the  form: 
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•  (2.i}-‘  (W(p))'r  "  rfr-bj)/  'n  TCr^Jldr 

c-i«  i“l 

i  i  7(nrP_1) 

•1,-lv,^  1  rn\ 


k 

I  I 


VgC«<-*  ;  ^ 


d-0  t  i^)6  C6(I)  C6(I) 

The  density  function  (6.3)  can  be  expressed  in  terms  of  Meijer's  G-function 
since  the  integral  in  (6.3)  is  expressible  in  terms  of  that  function  [10].  The 

density  function  of  W(p)  is  therefore: 


(6 


.4)  f(W(F))  -  frp(^n1+r*2))/rp(?n2)1 


k-0  k 


(t». ♦";))<  (?'1).  c«(!)  GP,0  m(p)  |*1 . *P) 

V*'”  p 


V  Y  Gr (W 

L  l  ^Ti  P.P 


1, 


l  l  2(nrp_1) 


d*0  6 

As  in  th.  previous  sections,  formula  (6  4)  -ill  (five  special  cases  as 


follows: 


a) 


If  we  let  0  -  0,  LY(0)  -  (in^  Cfi(I)  where  y  •  , 
and  after  making  use  of  (5.1)  we  obtain  formula  (4.7)  of  Pillai, 
Al-Ani  and  Jouris  [15]  for  testing  the  hypothesis  HQ:  yA  -  I. 

X  >  0  being  given. 


,«  Mu  a 
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b).  For  X  ■  1  and  A  ■  I,  and  after  making  use  of  (5.3)  we  let 

k  ■  6  *  o  (o  being  the  paratition  of  s  in  the  expansion  of  L^), 
then  we  have  formula  (5.2)  of  [15].  Note  that  in  this  caiie 
I  are  the  roots  of  determinantal  equation 

•  “  p 

^  -  ^(S^S^I  »  0,  and  is  the  Wilks'  criterion  for 

MANOVA. 

7.  The  density  function  of  the  largest  root.  In  this  section  we  derive  two 
expressions  for  the  density  function  of  the  largest  root  rp  of  S^1.  In 
obtaining  this  density  we  start  from  the  joint  density  of  the  roots  r^,...,  rp 
which  is  given  in  Pillai  [12]  by  the  formula: 

-tr  fi  -  hn  i(n.-p-l)  p 
CjCp.nj.nj)  e  ~|a|  ^Rr  II  (r^) 

1  1^1 
tr  Z  -  -sn.  7  i  T  "  7^ni4'no^ 

*  /  e  '|Z|  11  I  2~P| I+HRH' (I-Wj  A” 1 (I-W) Z |  1  1  dH  dZ 

Re  Z>0  ~  0(p)  .  ~  ' 

where  R,  Z  and  W  are  as  in  the  previous  sections, and 

Tp(p-l)  i  |p(P+l)  i 

C1(P»n1»n2)  »  1L  rp(T<Vn2))/[(2,ri)  y^)]  • 


Now  use  lemma  2  of  Khatri  [8]  by  taking  g(F)  ■»  rp  and  also  apply  formula 
(44)  of  Constantine  [1],  then  the  second  integral  in  the  above  expression 
becomes : 


12  1  11 > 
•5P  1  91  1  *  7^1 

[*  /rp(ip)]|i  ♦  (i-w)y1(i-wr  rpi 


.1  .1 

1F0(i(n1*n2);  r'^I  +  (I-W)  2A(I-W)  2  r"1)"1.  rp  I  -  R) 
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where  I  j  is  the  identity  matrix  of  order  p-1  and  I  is  that  of  order 


^p(p-i)  |p2  ,  4p(p+i)  t  , 

C2(p,n1(n2)  -  2  IT  rp(j(Vn2))/[(2,ri) 


Expand  jFq  in  terms  of  zonal  polynomials  and  for  integration  with  respect 

to  H  apply  lemma  3  of  Khatri  [8] ,  then  we  obtain  the  density  function  of  the 

largest  root  r  of  the  form: 

P 

-tr  n  -  •in  ^pn.-l 
(7.1)  C3(p,nltn2)  e  ~|a|  2  1  r2  1 


.  I  l  Sctla-jJ 

k-0  <  k!  rp  l(i(n1+P+l),K)  CJI) 


/  e  ~  | Z |  2  1|l  +  (I-W)  V1  (I-W) 2  r  |  2(W 

Re  Z>0  ~  . P 


•  CK[{I  ♦  (I-W)  2  A(I-W)  2  r” 1 } ~ 1  ]  dZ  , 


Let  hj^  ■  r^/Tp,  i  -  1,2,...,  p-1  and  H  ■  diag  (hj,...,  hpl),  then  the 

joint  density  of  r  ,  h. ,  h,,...,  h„  .,  where  0  <  r  <  •, 

P  l  2  p-l  p 

0<hi  <  h^  <  ...  <  h  j  <1  is  given  by: 


-tr  n  -  in  ipn  -1  i(n.-p-l) 

C2(p,ni,n2)  e  '|A|  rp  |H|2  |  Ip  j-H^II  (hi-hj) 


tr  Z  -  — ii  —  — 

/  e  '  |Z|  2  1 1 1  ♦  (I-W)  V1  (I-W) 2  r  I 
Re  Z>0  ~  . P 


1  11,, 

2,-1, T 


-  I  .  i 

’  iFo(I^nl+n2);  2  rp1)1*  ^-i"”)  dz  » 


PT 


f  w.fi 11 ,11W'I,UHWP»'I'  Ilflini"  1  I-B"! 


where 


C3(p,nl’n2)  "  C2(p,nl*n2)rp-l(i(nr1^rp-l(I(P"1)/n 


Now,  the  integral  in  (7.1)  can  be  written  as: 


L  ,>2 


tr  Z 


1 


.  -  4(ni+n->)  -  ,  -  i  -  4(n.*n-) 

/  e  ~|Z|  ?  1|(I-W)r  A_1|  2  1  2  |l+(I-W)  2r_1A(I-W)  2|  2  1  2 

Re  Z>0  '  "  '  P“  '--P*** 


C  ({I  ♦  (1-W)  7  r_1A(I-W)  2}-1]  dZ  . 

K  -  ~  ~  P  ~  ~ 


Applying  (12)  of  Constantine  [1],  the  above  integral  becomes: 


Ir.-.flW 

P.-_— _ 


rp(j(nl+n2).K)  Re  Z>0 


tr  Z  -  -  i(n  +n  ) 

/  e  ~  | Z  |  2  1 1 1-W|  2  1  2 


2  _-l.,T  2,,llc.l2(nl+VP"1) 


•/  exp[-tr{S(I*(I-W)  1  r’AA(I-«f)  ^)}]|s| 

s>0  ....  P  -  -  * 

•  Ck(S)  dS  dZ  . 

11  1, 

2  ~  2  i2^P+^ 

Let  B  »  (I-W)  S(I-W)  and  compute  its  Jacobian,  which  is  J(S;B)-|l-W| 
After  making  a  substitution  and  necessary  rearrangement,  the  above  expression  now 
can  be  written  as. 


1 


(7.2) 


|r  A 


-i  r  t'W 


4 


rp(j(nl+n2),  <)  n»0  v 


tr  Z  -  in. 

n!  L  6<  v  ;  e  ~IZI 

6  Re  Z>0 


11  ^  UtJ 


1, 


,  tr  B  vj-tn.+n ,-p-l) 

-  -  Ib4  1  Z 


/  exp[-tr{(I+r"  A)B}]  e  ~  | B | 
B>0  *  p  *  * 
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C  [(I-W)B]  dB  dZ  , 


i'  I  '  ,*i  lib 


where  y  are  constants  (see  Khatri  and  Pillai  [9])  and  I  6^  ■  k  +  n, 
6  ■  ).  The  second  integral  in  (7.2)  equals 


(7.3)  /  1 1 “  B|B|l(n.*vp-» 

B>0  a-0  a  4  L 


ca[(I-W)B]  dB  , 


where  now  (I+r  A)  is  a  diagonal  matrix.  Transform  B  +  H  B  I!'  where 
H  c  0(p),  then  the  volume  element 


(7.4) 


dB  -  n  (b.-b.)  n  db. (dH)  , 


where  b^  i  =  1,...,  p,  are  the  roots  of  B.  Making  substitution  in  (7.3) 
and  integrating  over  H,  (7.3)  now  becomes: 


,  ,  ;  r  (-l)a  '«<!<!>'.<!>  tr  !,.,I(nl*VP‘1) 

]■•■!  I  l  - irm -  e  l?l 


a=»0  a 


C.(I-K)  p 

•  ",  -c^ir  dbi 

i>i  J  6V,/  x»l 


where  now  B  is  a  diagonal  matrix  B  «»  diag  (b^...,  b  ), 


■  wsmm 


Now  let  us  transform  the  diagonal  matrix  B  back  to  the  symmetric  matrix  B 
and  use  (7.4)  then  integrate  out  H.  Then  we  have 

/  expf-trUl+r'^B}]  e^  ' jB|2(VVp‘i;) 

B>0  -  P  -  -  'J 

V»)  C6(I-W) 

C6(I)'  d?  • 

C  (I-W)  » 

'  Jo  l  *  l 

t  i(n1+n„-p-l) 

•/B>0  exP["tr{(I+rp1A)B}] J B | 2  1  2  Cy(B)  dB  , 

where  gJ>T  are  constants  and  Z  y.  =  Zfi.  +  t  ,  y  =  . 

Applying  (12)  of  Constantine  [1]  to  the  above  integral,  we  see  that  (7.2) 
becomes : 


.T  -1,'  I(ni+n2) 

I+r  A  2  n 

n~- -  y  y  tjjl  y  6  r 

(j(nl+n2).K)  n=0  v  n!  <5  K»v  ] 


L: _ E 

r 


Re  Z>0 


tr  Z  -  in 

®  ”l?l  C6(I-W)dZ 


'  1  Jo  l  J  J  <T  rp(Ttni*n2)-“>c„[(K1«'1l  • 


1  1 

Replacing  W  by  ft2  z"1  ft2 
combining  the  result  with  (7. 
root  rp  of  SjS'1 


and  making  use  of  (17)  of  Constantine  [2]  and 
1) ,  we  obtain  the  density  function  of  the  largest 


as: 


(7.5) 


-tr  ft  -  *n, _  ,  -  i(n,+n.)  ipn  -1 

C(p»i»1  ,i»2)e  “  | A|  1|l*rpA  l\  2  r2 


.  j  y  (|p»n.(l(p-n)K  jjta;,  .  3_ 

L  L  i  i  L  L  n  l  £  c  \) 

k*0  ic  k!  (y(n1^p^l))|c(«rp)K  n=0  v  6  #  (^ij) 


y(n  -P-1) 

l:  (n) 


6C«<P 


1^-1, 


■tI  I  tf  I  «5.,  <l‘V"2»p  cul(!*rp  ^  1  ' 

t=0  T  u 


where 


(7.6) 


C(p,n3  ,n2) 


r(T)rp(?(ni*r,2))rp-;(1t’*1) _ 

r (ip) r (Inj ) rp (in2) rp_  j (^(n  j  »P.i ) ) 


Note  that  in  obtaining  (7.5)  and  (7.6),  the  relations  r  (a,<)  =  (a)  T  (a) 

P  K  P 

and  C.(!p-l)/CK(V  ■  (y(P-l) )*/$?)*  have  been  used- 

The  density  in  (7.5)  will  give  (16)  of  Khatri  [8]  if  we  let  ft  =  0 
and  k  ®  6  ■  y. 

Alternately,  starting  with  (7.2)  and  directly  integrating  the  last 
integral  there  using  (12)  of  Constantine  [1],  we  get 


(7.7) 


- i - - -  l  l 

r  (T(ni+n2J»,c^  n“°  v 


tlL 

n! 


I  *<,v  rP4<V"2>-{> 


tr  Z  - 

*  /  e  -|Z| 

Re  Z>0 


1^ 

2  1  C6[(I-W)(rpA'1)]dZ  . 


Diagonalizing  (r^A  *)  by  an  orthogonal  transformation 
over  0(p)  and  finally  applying  (17)  of  Constantine  [2] 
function  of  the  largest  root  r^  of  S^S2*  as: 


H  and  integrating 
,  we  obtain  the  density 
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-tr  fl  -  icn.-l  * 


(7.8)  Ctp.n^n^e 


2*1 

rp 


k-0  <  kI  ('5'(n1+P'"1))<('51?)< 
i  -l  ^VP_1) 

y  y  (-1)"  y  »  )l» _ ^ 

-  "  “  -  ‘•v 


n-0  v 


where  C(p,nltn2)  is  as  in  (7.6). 


For  fi  ■  0  and  <  -6,  the  density  of  r is 


(7.9)  f(rp) 


2  w.  ,  r  > 

C(p,n  ,n_)  | A |  r  [  I - - T  *  ,1  “ 

1  2  '  P  k«0  k  k!  (y(n1*P+1))„;(2P)|C 


-  K  W1  * 


i  1  x 

-  C(p,n1,n2) |A|  2  1r|P  1  3F2(j(n1+n2) ,^p+l ,-(p-l) ; 

y(n1+P*l)»  ■jP!  rp^  ) 
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